Thermal noise for quantum state inference 
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In this work we describe a simple and efficient scheme for inference of photon number distribution 
by adding variable thermal noise to the signal. The inference remains feasible even if the scheme 
parameters are subject to random dynamical change. 
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Quantum state tomography is the most advanced and 
complete diagnostic tool available now. It allows one in- 
ferring maximal possible information about the state of 
physical system or about the process [ij. Quantum to- 
mography has already become a standard experimental 
tool enabling reconstruction even such fragile, exquisitely 
quantum objects as "Schrodinger cats" Intuitively, 
one expects that for performing measurements able to 
collect information sufficient to reconstruct a quantum 
state or process, it is necessary to build rather pre- 
cise measurement set-up. This intuition seems to be 
confirmed by existence of limits for detection efficiency 
required for performing a reconstruction, such as 50% 
threshold for a quantum homodyne tomography [3, 01 ■ 
Moreover, it seems natural supposing that it is necessary 
to know precisely what exactly one's measurement set-up 
detects (in more formal language, it is seems necessary 
to know all the elements of the positive valued operator 
measure, POVM, describing the experiment). 

Well, recently some works have appeared giving a 
rather obvious hint: with quantum measurements it 
might be really unnecessary to struggle for achieving ex- 
actly known (i.e. calibrated), and perfectly controllable 
measurement set-up. Under certain quite general condi- 
tions (such as, for example, Gaussianity of the state is 
question) it is possible to update information about the 
set-up and the signal state simultaneously (which was 
termed " self-calibration" ) [5j] . Very recently experimen- 
tal demonstrations of self-calibration were given Q. 
Moreover, it was demonstrated that sufficient knowledge 
about set-up can be acquired in the process of collecting 
data even in absence of any initial information about the 
measurement set-up [s^. 

However, self-calibrating approaches described above 
still suppose that there is some fixed measurement set-up, 
albeit possibly unknown one. But what if some noise is 
present? What if it is subjected to random (and possible 
uncontrollable) change? In this work we are demonstrat- 
ing that classical noise affecting the measurement set- 
up can be beneficial and usable for building robust, effi- 
cient and simple measurement set-up (which can even be 
much simpler than exiting methods) for certain tomog- 
raphy tasks. Moreover, even randomly changing POVM 
parameters might be not an obstacle for these tasks. 

A general impression about a possible role of noise 



can be given with the following simple example. Let 
us consider a set-up performing projection on the co- 
herent state with the amplitude a, i.e. with the POVM 
element 11 — |a)(a| (which can be realized with hetero- 
dyne measurement [9i]). If the amplitude a undergoes 
random changes (say, 5j) around some particular value, 
say, Q!Oj for a sufficiently large number of different 5j 
the resulting set of POVM elements (i.e. of the form 
IIj = |q;o + 5j){ao + 5j\) will be sufficient for performing a 
complete state/process reconstruction (which is attested 
by recent schemes of " data pattern tomography" 0, [lo| 
and quantum process tomography [Tl|). So, adding noise 
to a measurement can actually increase a region that the 
measurement set-up is actually "seeing", i.e. the search 
subspace. 

Now let us demonstrate how noise can be implemented 
for devising simple and efficient set-up for inference of 
photon number distribution of a single-mode state of elec- 
tromagnetic field. Notice that, for example, for quantum 
homodyne tomography an inference of photon number 
distribution is not much simpler than inferring the com- 
plete density matrix (one needs performing the same set 
of quadrature measurements). The task can be made 
easier by building the specific set of POVM elements. 
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where |n) are Pock states with n photons. In the end of 
90-s it was demonstrated that implementing only a single 
bucket detector with a set of variable absorbers changing 
the efficiency of the detection it is possible to collect data 
sufficient for inferring the photon number distribution 
p^ . For this set-up one has r„ = (1 — 77)", where 77 is 
the efficiency of the detector. The method was realized 
experimentally [l^. It was shown also that adding a 
coherent shift to the signal state it is possible to perform 
a complete tomography p^ . However, an experimental 
realization of the scheme remained rather challenging due 
to necessity to perform calibration of absorbers for signals 
of a few-photon level. To avoid this more sophisticated 
set-up was suggested and realized; there the signal travels 
along the fiber loop splitting on each pass, seefTEI. [l6j). 

Implementing noise gives a way to avoid using both 
sets of calibrated absorbers or loop detection for the task. 
Indeed, let us simply add a thermal noise (i.e. a source 
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FIG. 1: An example of two-photon Fock state maximal like- 
lihood estimation using variable thermal noise. 30 values of 
average number, n, of thermal photons were used equidistantly 
distributed in the interval [0.1,0.95]; the search space was from 
to 12 photons, 10,000 iterations of the reconstruction algo- 
rithm were taken, starting with the maximal entropy state. 
The detector efficiency, rj, is 0.8. The inset shows a depen- 
dence of fidelity on the number of iterations, L, for 20,000 
measurements (solid line) and 10,000 measurements (dashed 
line) per each value of n 



of detector "dark counts") with the average number of 
photons n to our signal. As follows from the Mandel's 
formula, without the thermal noise mixed with the signal 
the probability to have no clicks on our bucket detector 
is given by 



piv) 
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where /9„„ are diagonal elements of the signal state den- 
sity matrix in the Fock state basis; notice that p{ri) is, in 
fact, the photon-number generating function. When the 
thermal light with the average number of photons, n, is 
added to the signal at the entrance of the detector, the 
probabihty ^ is modified as to [Tt I fli ] 
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As it follows from Eq.(|31), one can just mix thermal 
noises corresponding to different temperatures to build 
the POVM elements set required for inference of photon- 
number distribution. In Fig|T] an example of the two- 
photon Fock state reconstruction is given. For the task 
30 different values of the average number of photons of 
the added thermal noise are taken. Efficient expectation- 
maximization iterative algorithm for maximal likelihood 
estimation was implemented [l^ [ij]. For 10,000 mea- 
surements per each value of n a fidelity of more than 
95% was reached for 10,000 iterations starting from the 
maximal entropy state. One should point out two im- 
portant features of the suggested scheme. The first one 
is necessity to use rather large search subspace due to 
the photon-number distribution of the thermal state be- 
ing long-tailed. The second one is the slower convergence 
of the algorithm for larger number of measurements (see 



inset in Fig[T]). Increasing of the number of measure- 
ment ultimately is eventually leading to increase in the 
reconstruction accuracy. However, with larger number of 
measurements one can actually get worse results for the 
same number of iterations. 

It should be stress out that adding thermal noise leads 
to significant simplification of the process of building the 
POVM elements set required for the reconstruction. In- 
deed, one is not even obliged to actually change param- 
eters of the source of thermal noise; it is sufficient to 
change a time-window of detection to change effectively 
an average number of thermal photons. Notice, that ther- 
mal noise needs not to be pre-calibrated: it is possible to 
change noise temperature arbitrarily and calibrate it by 
temporarily switching off the signal and collecting data. 
It means that one can actually perform the reconstruc- 
tion with usual daylight. The only obstacle is dispersion 
of the detection efficiency, so it is necessary to filter ther- 
mal noise to provide for the constant detector efficiency 
in the spectral interval of thermal noise and the signal. 

Of course, when one tries to build set of POVM by 
adding noise, the question arises of its influence on pos- 
sible reconstruction errors. Generally, the problem of er- 
ror estimation for quantum tomography is rather "hot" 
and controversial subject nowadays. Both method bases 
on ML estimation and Bayesian inference were suggested 
for the purpose (see, for example, Ref . |20l] and references 
therein). However, for diagonal elements inference with 
general POVM ([T]) we can suggest a simple estimation 
of error of the ML method along the lines suggested in 
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(and somewhat in the spirit of approach used in 



Let us consider our measurement with the general 
POVM ([T|) as the set of measurements with K complete 
POVMs each having just two elements Hj and / — H 
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j = 1, . . . , /-iT . In the limit of large number Nj ^ 1 of runs 
with the j-th POVM of the set, the likelihood function 
has the following form 
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with the number of outcomes Sj corresponding to the 
POVM element Hj (and, respectively, Nj — Sj for I—Hj), 
can be approximated by using the well-known large- 
limit of the binomial distribution (see, for instance, Ref. 
[2l|). The latter is a Gaussian approximation, which in 
our case gives 
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with 



a] = tr([/ - nj]p)tr(H,p). 
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then for the probability one has simply 



FIG. 2: An example of two-photon Fock state maximal like- 
lihood estimation using random variable thermal noise (a). 
The noise was taken to be Gaussian; 30 values of average 
number, n, of thermal photons were used equidistantly dis- 
tributed in the interval [0.1,0.95]; the variance is 0.1; the 
search space was from to 12 photons, 1,000,000 iterations 
of the reconstruction algorithm were taken, starting with the 
maximal entropy state. The detector efficiency, rj, is 0.8; the 
number of measurements is 50,000 per each value of n. The 
inset shows a dependence of fidelity on the number of iter- 
ations, K, for 50,000 measurements and variance 0.1 (dash- 
dotted line); 10,000 measurements and variance 0.01 (dashed 
line); 10,000 measurements and variance 0.1 (dotted line) per 
each value of n. The panel (b) shows an estimation of max- 
imal errors made according to Eq.® for different number of 
trials, A'^ x 1000 per average value (solid line with dot marks) , 
in comparison with the estimation for non-random POVM for 
fixed noise values equal to average ones. For this figure the 
average number of thermal photons is n = 0.1; the variance is 
0.1; the search space was from to 12 photons (dotted line). 



For sufficiently large Nj the width of the peak around the 
maximum likelihood point, i.e. the admissible values for 
Pmm, m = 0,1, . . . , M, is given by the variance appearing 
in Eq. ([5]). Therefore, the error S for each pmm of the 
maximum likelihood estimate is on the order 
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An example of maximal error estimation for the two- 
photon Fock state reconstruction is shown in Fig. [IJb). 
One can see that Eq.([6|) gives rather reasonable estimate 
of a maximal error for different number of measurements. 

The most obvious prerequisite for the "thermal noise" 
reconstruction discussed above seems to be a constant 
predefined level of the noise. However, it is not hard to 
see that noise can be varied during the data collecting 
process. It can be even random. Indeed, if the set of 
parameters, fl, describing a particular POVM element in 
Eq.©, n(^), represents a continuous random variable. 



(7) 



Of course, for each run of the experiment the result will 
depend on the particular value of random parameters, 
fl. The reconstruction procedure hinges on the fact that 
for the sufficiently long sequence of trials one can assume 
the all results were obtained with the averaged POVM 
elements ([7]). For our example of the "thermal noise" 
reconstruction these average POVM elements are to be 
defined on the calibration stage preceding a measurement 
of the state mixed with noise. Figl2][a) demonstrates 
that the reconstruction with "noisy" POVM elements is 
indeed feasible. There an example of two-photon state 
inference is shown for the "thermal noise" reconstruc- 
tion scheme with average number of photons fluctuating 
with the normal distribution. It is remarkable that one is 
able to achieve quite accurate reconstruction results with 
rather strong noise (when the variance of noise is com- 
parable with the average value of n) . The pay-off is the 
necessity to increase the number of measurements per as- 
sumed averaged POVM element. However, this increase 
is not crucial (for instance, just five times for the example 
shown in FiglJl corresponding to quite large variance of 
the noise). An estimation ^ shows that even for mod- 
erate number of trials maximal errors for fixed thermal 
noise and random one can be quite close (this situation 
is illustrated in the Figl2l[b)). Also, as the inset in Figl2] 
demonstrates, random variation of the thermal noise do 
not noticeably worsen convergence of the reconstruction 
procedure in comparison with the "fixed noise" results 
shown in FigHJ 

So, we have established that random variations around 
some fixed values of POVM parameters do not spoil the 
reconstruction. Now let us demonstrate that we can re- 
lax our requirement for controlling the experiment even 
further. Values of POVM parameters might be not fixed 
at all. A knowledge of average values of these parame- 
ters in a certain time interval is sufficient for performing 
the reconstruction. In Fig|3] an example of the recon- 
struction is shown for the average number of photons of 
the mixed thermal state undergoing the stochastic drift 
process. Since copies of the signal are assumed to be 
generated as an equidistant sequence of pulses, random 
change of the average number of photons was taken to 
be described as ID random walk process. Probability of 
jump to higher n was taken to be slightly larger, than to 
lower n (1% higher), so the average number of photons 
was gradually increasing (see FiglS]). Consequent sets 
of 30,000 measurements were taken for averaging and as- 
signing an averaged POVM element (whole of 30 different 
values of n non-equdistantly distributed in the interval 
[0.1,0.952]). As it can be seen in FiglSj this approxima- 
tion has allowed to achieve rather good estimation of the 
two-photon signal state even with rather moderate num- 
ber of measurements. 
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FIG. 3: An example of a realization of the randomly chang- 
ing average number of thermal photons. Random change is 
ID random walk with probability 0.51 of jumping to a higher 
value of n (they are taken to be a discrete set with the step 
5 * 10"^. The total of 9 * 10* measurements were assumed. 
For the purpose of the reconstruction this set was divided on 
30 subsets with 30,000 measurements in each. So, 30 differ- 
ent POVM elements ^ were assumed for the reconstruction 
with different n corresponding to averages over each subset. 
10® steps of the iteration procedure were takes. The inset 
shown the convergent result of the estimation of the signal 
two photon state.) 

Concluding: we have established that noising the mea- 
surement (or adding noise to the signal) can be a simple 
and efficient way to produce a set of POVM elements 
sufficient for reconstructing the state of the signal. Mix- 



ing the signal with the thermal noise on the detector (i.e. 
adding "dark counts"), one can make very simple and 
non-expensive set-up for estimating a photon number dis- 
tribution. One does not require for it a set of calibrated 
absorbers or loop detectors. Calibratioir of the thermal 
noise can be done directly in the process of measurement. 
Randomness of the added noise does not spoil the recon- 
struction, provided that the average values of noised pa- 
rameters are known for sufficiently long series of trials. 
We have demonstrated that both for random parameters 
fluctuations around some set of fixed values, and for the 
random change similar to stochastic drift/random walk. 
The latter feature hints at the possibility to overcome 
the main reason of using the same source for generating 
the supposedly unknown signal and the reference field in 
the schemes of complete quantum tomography: the phase 
drift. Thus, this work is the step to reaching the ultimate 
goal of a quantum state tomography: accomplishing the 
reconstruction of truly unknown signal state. 
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